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SUMMARY 

We consider the effects of critical layer nonlinearity on spatially grow- 
ing oblique instability waves on nominally two-dimensional shear layers between 
parallel streams. The analysis shows that three-dimensional effects cause non- 
linearity to occur at much smaller amplitudes than it does in two-dimensional 
flows. The nonlinear instability wave amplitude is determined by an integro- 
differential equation with cubic type nonlinearity. The numerical solutions to 
this equation are worked out and discussed in some detail. We show that they 
always end in a singularity at a finite downstream distance. 


1 . INTRODUCTION 

External excitation of (convecti vely unstable) free shear layers between 
parallel streams produces spatially growing instability waves that are ini- 
tially governed by linear dynamics for sufficiently small excitation ampli- 
tudes. While the instability amplitude continues to increase with increasing 
downstream distance, its local growth rate must ultimately decrease due to vis- 
cous spreading of the mean shear layer. Nonlinear effects can then become 
important in a "critical layer" at the transverse position where the mean flow 
and instability wave phase velocities are equal (once the instability wave 
amplitude becomes sufficiently large and its growth rate becomes sufficiently 
small). The unsteady critical layer flow is then governed by a nonlinear 
vorticity equation, while the motion outside the critical layer remains essen- 
tially linear. The external instability wave growth rate is, however, com- 
pletely controlled by the nonlinear dynamics of the critical layer. 

There are now too many nonlinear critical layer analyses in the litera- 
ture for us to summarize here. We refer the reader to the excellent review by 
Maslowe (1986) for nonlinear critical layers in general and to the one by 
Stewartson (1981) for Rossby-wave critical layers in particular. Here it is 
only appropriate to mention the analyses specifically concerned with spatially 
evolving flows. Huerre (1980, 1987) considers the two-dimensional incompres- 
sible shear layer in the viscous critical layer regime where the nonlinear 
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effects enter only at higher order. Huerre and Scott (1980) and Robinson 
(1974) consider the strongly nonlinear regime, but their choice of scaling pre- 
cludes the possibility of matching their results onto the linear upstream 
solution. 

Goldstein and Leib (1988) and Goldstein and Hultgren (1988) overcame this 
difficulty, but again consider only the incompressible case. They consider 
only a two-dimensional flow, since the two-dimensional instability wave is the 
most rapidly growing linear mode in that case. However, the oblique modes 
exhibit the most rapid growth in sufficiently high Mach-number supersonic shear 
layers (Gropengei sser , 1969; Jackson and Grosch, 1988). Goldstein and Leib 
(1989) consider the case where the unsteady flow evolves from a single oblique 
mode, in which case it is appropriate to suppose that the initial instability 
wave grows in its propagation direction. The cross flow velocity effects can 
then be eliminated from the analysis by use of an appropriate "Squire trans- 
form," but the critical layer nonlinearity still behaves quite differently 
from the incompressible case - primarily because the temperature fluctuations 
have algebraic singularities in the critical layer and therefore become very 
large relative to the remaining velocity components causing the critical layer 
nonlinearity to occur at a much smaller amplitude vis-a-vis the two-dimensional 
isothermal case. The resulting critical layer nonlinearity is then weak in the 
sense that the flow is governed by linear dynamics to lowest order of approxi- 
mation, with nonlinearity entering only through the higher order (inhomoge- 
neous) terms. The instability wave growth rate is still controlled by the 
nonlinear terms, but can now be calculated from an amplitude equation similar 
to the one found by Hickernell (1984) for the Rossby wave singular modes. 

Here, we consider the case where the initial instability wave grows in 
the streamwise direction. It is then appropriate to suppose that there are 
two oblique modes with the same frequency and streamwise wave number but with 
equal and opposite (real) spanwise wave numbers in order to represent a fixed 
spanwise structure. (In real flows the allowable spanwise wave numbers might 
be selected by the side-wall positions.) The cross flow velocity fluctuations, 
which have the same algebraic singularity in the critical layer as the tempera- 
ture fluctuations, now become coupled to the velocity fluctuations in the plane 
of the wave, causing the critical layer nonlinearity to again become important 
at smaller amplitudes than in the two-dimensional incompressible case. In fact 
the nonlinear oblique mode interaction causes the critical layer nonlinearity 
to occur at even smaller amplitudes than in the single mode compressible case 
(Goldstein and Leib, 1989). Critical layer nonlinearity now becomes important 
when the instability wave growth rate is 0(e^3) rather than 0 , as in 

the Goldstein and Leib (1989) case, where e denotes the characteristic insta- 
bility wave amplitude at the start of the nonlinear region. 

While the phenomenon is of most importance in supersonic shear layers, we 
restrict the analysis to the incompressible case because the nonlinear criti- 
cal layer dynamics are unaffected by compressibility effects, and their inclu- 
sion would greatly complicate the analysis. Moreover, the final amplitude 
equation (3.69), which is the principal result of this paper, remains valid 
for the compressible case. The instability amplitude outside the critical 
layer can again be determined from an amplitude equation because the critical 
layer nonlinearity is still weak. The equation is similar to that of Goldstein 
and Leib (1989) in that it is an i ntegro-di fferenti al equation with cubic non- 
linearity, but the structure of the nonlinear kernel function is now somewhat 
different from theirs. The equation still has to be solved numerically. This 
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is accomplished by using a fourth-order predictor-corrector scheme to inte- 
grate in the downstream direction, starting from the upstream linear state 
which is prescribed far upstream in the flow (relative to the streamwise 
length scale of the nonlinear region). 

The calculated instability wave amplitudes initially follow the pre- 
scribed linear growth, but soon begin to either saturate or increase their 
rate of growth when the nonlinear effects come into play. Cumulative history 
effects eventually reverse these trends, causing a rapid increase in amplitude 
which ends in a singularity at a finite downstream distance. The local asymp- 
totic solution of the amplitude equation is derived. It. shows that the flow 
will become fully nonlinear everywhere in the shear layer and that the motion 
will then be governed by the full three-dimensional Euler’s equations in the 
next stage of evolution downstream of the weakly nonlinear region. 

The overall plan of the paper is as follows. The problem is formulated 
in section 2, where we show how the nonlinear critical layer gradually evolves 
from the strictly linear finite growth rate solution. The flow outside the 
critical layer is a linear, unsteady three-dimensional perturbation about the 
two-dimensional mean shear layer flow, which can be treated as locally paral- 
lel on the streamwise length scale over which the nonlinear effects take 
place. The latter occur entirely within the critical layer to the order of 
approximation of the analysis and determine the unknown amplitude function in 
the external solution. The transverse velocity fluctuation is taken as the 
basic variable for the external flow, and the (linear) equation for this quan- 
tity is expressed in terms of the "Squire coordinates" associated with either 
of the two oblique modes. 

We introduce a "slowly varying" amplitude function, which depends on the 
streamwise coordinate through an appropriately scaled variable and is ulti- 
mately determined by the nonlinear flow in the critical layer, which we analyze 
in section 3. Matching with the linear external instability waves leads to the 
nonlinear integro-differential equation that determines the amplitude of those 
waves. The numerical and asymptotic solutions to this equation are discussed 
in section 5. 


2. FORMULATION AND SOLUTION OUTSIDE THE CRITICAL LAYER 


We are concerned with a nearly inviscid incompressible shear flow between 
two parallel streams with nominally uniform velocities U (1) > U ( 2>. The 
upstream flow consists of a steady two-dimensional shear layer and a pair of 
oblique (i.e., three-dimensional) spatially growing (i.e., time harmonic) 
instability waves with the same frequency and streamwise wave numbers but with 
equal and opposite (real) spanwise wave numbers. The streamwise, transverse, 
and spanwise coordinates (x, y, and z, respectively), the time t and all 
velocities are normalized by 6 0 , S 0 /A, and A, respectively, where 6 0 is 
half the mean shear layer momentum thickness and 


is a measure of the velocity difference across the shear layer. When nonlinear 
effects do not first intervene, the gradual viscous spreading of the mean shear 
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layer causes the common spatial growth rate of the two linear instability waves 
to gradually decrease until they approach their common neutral stability point, 
whose Strouhal number and streamwise and (real) spanwise wave numbers we denote 
by S 0 , a, and ±3, respectively. 

As already indicated, nonlinear effects first become important at the 
streamwise position (upstream of the linear neutral stability point) where the 
local Strouhal number S (frequency normalized by A/S 0 ) is 

S = S Q + e ^ /3 S j . (2.2) 

Here S] < 0 is assumed to be 0(1) constant, and e denotes the characteris- 
tic amplitude of the instability wave in this region. The instability wave 
growth rate, which is also 0(e'^), will then be determined by the nonlinear 
critical layer effects. Finally, we require that the origin of the spatial 
coordinates x, y, z be located within this nonlinear region. 

The transverse velocity fluctuations of the two upstream linear instabil- 
ity waves will then be of the form 



(aC-e 1/3 S ] t)-l/2 S ] U^kX 1 


where t denotes the time, 


a 



(2.3) 


±3 denote the (real) transverse wave numbers of each of the individual insta- 
bility waves and ±0, where 


9 = tan 1 — 


a 


(2.4) 


denotes the direction of propagation (relative to the mean flow direction) of 
these waves at their common neutral state. 


5 E X - U c t 


(2.5) 


is a streamwise coordinate in a reference frame moving downstream with the 
neutral phase velocity 


U £ 3 


x ] - C 1/3 X 


( 2 . 6 ) 

(2.7) 


is a scaled streamwise coordinate (in the fixed reference frame) and, 

i s ii • — 

Vc 2 


is a scaled complex wave number whose imaginary part is minus the common growth 
rate of the linear near neutral instability waves. <t>j^(y) is an appropriate 
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solution of Rayleigh's equation, which can be taken as the neutral solution to 
the required order of approximation. The complex constant a> is a measure 
of the complex (scaled) amplitude of the two waves. 

The normalized pressure is denoted by p and the normalized velocity 
components in the moving {£, y, z} coordinate system by {u, v, w}. We expect 
the solution outside the critical layer to expand like 

U = U(y) + eu ] + e 4/3 U 2 + . . . (2.8) 

v = ev ] + e 4/3 v 2 + . . . (2.9) 

4/1 

w = ew^ + e w 2 + . . . (2.10) 

p = ep.| + e 4 ^ 3 p 2 + . . . (2.11) 

where U(y) + U c is the base-flow velocity and the form of the upstream linear 
solution suggests that V] will be of the form 

v i = -a(e 1,3z + e~ 1|32 )& i4>{ 1) (y)A t (x 1 )e la ^ , (2.12) 

where 


C = K - S lC 1/3 t/a , 


(2.13) 


and At is a function of the slow streamwise variable x i , which will ulti- 
mately be determined by the nonlinear flow in the critical layer, but matching 
with the upstream linear instability waves requires that it satisfy the 
upstream boundary condition 


a^-S.U^x 


/ 2 


as 


(2.14) 


The 0(e 4/ 3) terms (in (2.8) to (2.11)) are induced by nonlinear effects within 
the critical layer and by derivatives of the 0(e) terms with respect to xj. 
They are, therefore, at least in part, associated with linear growth rate 
effects. In fact, the entire solution (2.9) satisfies linear dynamics to the 
indicated order and is therefore determined by Rayleigh's equation 


(ft * u kV ' 1 - u " I - 0 


(2.15) 


where the primes denote differentiation with respect to y and 


V 





(2.16) 


is the Laplacian in the moving reference frame. It follows that V 2 is a 
function of the form V 2 (C, *1, y, z) and therefore satisfies 
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(2.17) 


(UV 2 




8v l [fl 3_ 

9x^ + a 9£ 


(U 


V 


_9 

9x 


V 2 v , 


where now denotes the Laplacian with respect to £, y, and z and the 
partial derivatives with respect to t; are at constant X] and vice-versa. 
It follows from continuity that 


9u 2 9v 2 9w 2 9u 1 

9C~ + 9y~" + 9z~ + 9x~^ = 0 


(2.18) 


4> 


(1) 

1 


must, of course, satisfy the reduced Rayleigh's equation 


L <1> *j 1> = 0 (2.19) 

where 

L<'> i t |f) <2.20, 

dy z v u J 

denotes the reduced Rayleigh's operator, and and w^ are then given by 

u | = (e 1,5z + e _1|3z )cos e^«[*j 1> '(y) + tan 2 e ) (y)]A t (x ] )e 10 ^ 

+ ^F^ ) (y,x 1 )e 2l ^ z , (2.21) 

and 

w 1 = -2 sin 0 sin 4>j n - <J > (1 ) ‘)aV°^ , (2.22) 

where the C-i ndependent term in (2.21), which is governed by linear dynamics 
in the main shear layer, is induced by nonlinear effects in the critical layer. 

Then since U 2 , V 2 , ancl w 2 have the same £-periodi c i ty as V], they 
must be expressible in the form 

oo 

u 2= ^X F 2 m) < — i )e ’^ (2.23) 

m=0 


v 2 = -a Re i ^ $ 2 m) (y,z,x 1 )e imai; (2.24) 

m=0 


w 2= ^X H r<y-i>e^ <2.25) 

m=0 
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Substituting (2.12) and (2.24) into (2.17) we find that 



= 2 


cos |3Z <t> 


(1) 

2 


(y.Xj) , 


(2.26) 


where <t> 


(1) 


sati sf ies 


laUL"^” 


iS l - U c fx 




(AW -2 A <1>\ , 2 ll dA t .(l) 0 ,,, 

1*1 " a *l + 2ot U dlTT ^1 (2 - 27) 


It now follows from (2.23) to (2.25) and (2.18) that 


.( 1 ) 


(sec + sin 6 tan 9 H^^e 1 


i*z + e -1»z) 


(2.28) 


and 


H< n . -sin sh' 0 ^ 132 - r IBz ) 


(2.29) 


where F- 1 ^ is given by 


r ( 1 ) ,(1)' cos 2 e {.<])• . 2 ft U' (lA dA 1 

F 2 = *2 --UT" *1 +tane 0"*l | d)T 


(2.30) 


and H 2 can be defined similarly, but will not be needed in the following 
analysis . 

Equations (2.19) and (2.27) must, in general, be solved numerically. 
Fortunately, we need only know the local behavior of their solutions at the 
critical level where U(y) = 0. This will occur at the inflection point 
U"(y c ) = 0, or 

U“ = 0 , (2.31) 


(where the subscript c is used to denote quantities at the critical level) 
provided we assume, as we now do, that U(y) has only a single inflection 
point. 

The critical level, which we can always suppose to lie at y = 0, is then 
a regular point for the operator (2.20) and equation (2.19) will, therefore, 
possess the two linearly independent solutions 


•r(i) 

<f> 



(2.32) 
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and 


$ 


( 2 ) 


y + 


1_ 

3! 



(2.33) 


as y - 0. 

Then we can put 


<D 


( 1 ) 


r<l> 

= $ + 


h * (2) 

b^<p 


(2.34) 


where b] is a constant which must, in general, be determined along with a 
by solving (2.19) numerically subject to the proper boundary conditions at 
y = ±oo. Then (2.21) and (2.22) imply that 


u 1 = cos 9 


( el 


i (3z „-i(3z\ j tan 2 9 


+ e 


^ 2 ( 


+ b^l + tan 9 


LI' ' ' 


1 ”c (, 5 . _2 


1 + 2 tan 9 




1 + | tan 9 


y 


+ /?eF;° ) (0,x 1 )e 2i3z 


+ . 


(2.35) 


w^ = -2 sin 9 sin |3z 


V] = -a(e i|5z + e- if 3 z )^iA t e i<xC ♦ . . . 

&MA + e iaC + . 


i i u i 

y " 2 a + 3 Ul 


(2.36) 


(2.37) 


as y > 0, where we have anticipated the fact that the critical layer solution 
produces an 0(e) mean flow change across the critical layer but no change in 
the second harmonic of the streamwise velocity fluctuation. 


It follows from (2.27), (2.32), and (2.33) that there exist two continuous 
functions, say 4>pj and <J>p 2> which satisfy 


, ( 1 ). U" ( 1 ) 

L * p j - -J ^ 

’ 1 U 


(2.38) 


and 


L (1) <t> - t> (1) 

L P,2 1 


(2.39) 


but are, in general, unbounded at y = ±®. T>pj will behave like 

ei + e2,Ly lyl + ®2y + • • • , as y -* 0 and <t>p2 will be regular there. 

The relevant solution to (2.27) must then be of the’form 
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(2.40) 


^’■1 K - ,S 1 A+ X* P.l - C 2.1 ?<,> - b 2,l J<2 ’. 


- 2 iaA^ ( <t> - c* 4 >^ ^ _ 5 * 

^ 1 ^' X 1V P,Z 2,2 2,2 


where bf n* c 2 n are real constants (even on the slow scale x i ) , which must, 
in general, be determined numerically. 

It now follows from (2.30) that 



sin 2 9 .t 
lay X] 


+ e. 


+ e ? In |y| + 




- a ( U C A I, - 1S 1 A+ J b 2.1 * ' ' ' JS J 1 -* 0 (2 ' 4,) 

where the specific forms of e] and e 2 are immaterial, and we do not need 
the corresponding expressions for in the following analysis. 


3. THE CRITICAL LAYER 

Equations (2.35) and (2.37) clearly show that the outer expansion (2.8) 
to (2.11) becomes singular at the critical level. The linear small growth 
rate critical layer thickness is of the order of the linear growth rate which 
is 0(e 1 ' 3 ) in the present case. The appropriate scaled transverse coordinate 
in this region is therefore 


Y = 


1/3 


(3.1) 


Introducing this along with (2.12), (2.23) to (2.25), (2.40), (2.41), and 
(2.32) to (2.37) into (2.8) to (2.11) and reexpanding the result suggests that 
the critical layer solution should expand like 


1/3 2/3~ ~ 4/3^ 

U = c / u ' Y + e / u A + eU, + e ,J U, + . . . 
C 0 1 2 


V = -C 


2/3- 


a(e i(5z + e-’^&lAV* ♦ ♦ « 4/3 I 


2 + 


2/3~ 


4/3' 


W = e Wq + eW-j + e + . . . 


4 / 3 : 


5 / 3 : 


p = cp Q + g P 1 + e P 2 + 


(3.2) 

(3.3) 

(3.4) 

(3.5) 
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where we suppose that the e n Sin e terms have been Incorporated in the uq, 
u 1 , etc . , 

S 0 . U' cos e(e i8z , e' 1 ® 2 ) &*♦.'“< , (3.6) 

and we have put 

v = e ~ 1/3 v (3.7) 

for convenience. 

The u n . v n , w n , etc. are functions of C, Y, z, and X], only and are 
determined by the inviscid momentum and continuity equations. It follows from 
(2.23), (2.25), (2.28), and (2.41) that matching with the external solution 
requires that 


P 2ir/a 



f 2ir/(3 

e -l(aC+Bz ) AQ +^ >z>x ^ d2 dc 

J 0 



- -fu A f 

a\ c x. 


- iSjA^b 


2,1 



(3.8) 


where we have put 


and 


Au~ - 1 im 

^ Y+co 


<C»Y,z,x 1 ) - u^(C ,-Y,z, ) 


u„ = u„ cos 9 ± w sin 0 
n n n 


for n = 0,1,2, 


(3.9) 


(3.10) 


The requirement (3.8) merely states that the change in propagation 
direction velocity across the critical layer, as calculated from within, is 
equal to the change in that velocity as calculated from the liner external 
sol ution. 


It turns out to be convenient to work in terms of the spanwise vorticity 


_ 8v _ 8u 
8x 3y ’ 


which expands 1 ike 


“ = -u’ 


1/3; 2/3~ 

e Uqy - e U-| Y + + . . . , 


(3.11) 


(3.12) 
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where 


&>2 = ~ u 2y + aa ' e + e 


( e is z + 


(3.13) 


Then to and 


satisfy 


„ 1/3 11 3_ - 3_ 

U * e a /3C + 3Y + e 


U = {u.v.w} 


1/3 


( U + U c) dt + w h 




1 


1/3 1 i3_ - 3 1/3 

U-c — ^7+V^n+e 

a jaL, 3Y 


(u . U c ) 


3 u 3_ 
8x 1 + 3z 


u = 


( l/3„ _JL „ „ I 

J P C + e c 2/3 P Y* p zj 


and 


u ? * 5 Y ♦ w 2 * e l/3 u . 0 


(3.14) 


to * -jft u z w Y ♦ «w z 


(3.15) 


(3.16) 


(3.17) 


Substituting the expansions (3.2) to (3.6) and (3.12) into (3.15) and (3.16), 
we find that 


W Q = -u^a sin e cos 9 Qe\(e^ z - e 1 ^A^e 1 

V, . <.•* , ,-"*)(& ,aV«% - & . s 0 fe ); 0 - u.y; 0X) - ;, 2 


(3.18) 


^ 0 u Q = sin 2 e(e 1Pz ♦ e'^/We 10 * 


(3.19) 


(3.20) 


Vir - <° ,BZ * *-<•“)( w^)c 0YY . (v 02 - = 02 ; 0 ) y ♦ u '(; l2 - % x Y ) 


(3.21) 
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*0*2 


and 


S (e' B2 ♦ e- iBz )(/? e M + e '^)(5; Yy - cos «”■») - -§T-(e l(!Z ♦ e' iB2 ) 
X /?«(u c aJ - lS,A^)e la ^ ♦ (g-; 0z - u 0 u| c - ; 0 g; 2 - cos eu 02 


- “l sl " e “oc - *l“ov) y + U c( cos 

- (v*ox,) Y 


0w o _ - sin 6w or - sin 0w^ - Yu 


2z 


2C 


j Ix,y) 


i '"i 
( 3 . 22 ) 


V 1Y " U K W lz ■ u 0x, ’ 


(3.23) 


where we have put 


- U -d— + [ u 1 Y - — |— 
■*0 - U c 9x 1 y U c T a J9C ’ 


and 


q, E -w, 


cos 0 + sin ©(^Y + 2a|3 sin |3z Re iA^e’ a ^ 


^(e ,Bz ♦ e- 1 Bz )/?«A t e ,a? 


(3.24) 


and matching with the outer solution shows that q£ and u^y - 1/2 U^.' 1 cos 0Y^ 
go to zero as Y -* +«> . 


It is now convenient to introduce the following normalized variables 

x - “ 2 S l U c x l ” x 0 ’ 


(3.25) 


n = -2al Y aU , / s i u c * 


X E a( - X Q , 


(3.26) 

(3.27) 
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and 


A = 4A t a 2 e 0 ( Sl U c ) 2 U- , 

where xq and Xq are real constants. 

Then (3.18) and (3.20) can be Integrated immediately to obtain 

u 0 = - - C - 1 sin 0 tan 0 cos |3z #0iW^ O) e 1X , 
u a l 

- Mr S 1 , v 

w Q = - q ~ sin 0 sin |3z Re wj 0) e lX , 

where we have put 

W 1 (0> = J X e _1n(x_x> A(x)dx . 

— 00 


It is worth noticing that 



(3.28) 


(3.29) 

(3.30) 


(3.31) 


when 


A -» e^x 

so that (3.29) approaches the linear critical solution (given by the generali- 
zation of eq. (4.34) in Goldstein and Leib 1988) when AT approaches the lin- 
ear upstream condition (2.14). 

Equations (3.19) and (3.21) clearly possess solutions of the form 

2 

u ]Y = \ U“'Y 2 - \ s1n eu c ^ Q n 1 | !i (n ’ x) e 1(nX+mBZ) » <3.32) 

n , m=— 2 

and 

2 

*i = ^ S lUc U' sin 0 T W^ 1 J i ( n ,x)e l(nX+m,3z) , (3.33) 

n 

where the first term in (3.32) was inserted to insure that ^ 0 as 

n +® for all n,m. n * m 

Substituting (3.32) and (3.33) along with (3.29) and (3.30) into (3.19) 
and (3.21) shows that 
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lX” = i sec e(2 Sin 2 e 


,(0) 


2 - jUA*w{° ) ) 


I , 1 > A 

L 2 W 2,0 " 0 


L 2 Q<!> = tan e sec 9^ - 41 sl^ew' 0 '^’) 


(3.34) 

(3.35) 

(3.36) 

(3.37) 


where the star denotes the complex conjugate, we have put 

L E -r-+ inn for n = 0,1,2, . . . , 
n 9x 

and in writing down (3.36) we have anticipated the fact that the relevant 
solution to (3.35) is the trivial solution 

= 0 . (3.38) 

Inserting (3.31) into (3.34) and (3.36) and integrating yield 

,x 


Wq 1 ^ = i sec 0 Re 


A*(x)e" 1nX [l 1 (x) + 4 si n 2 0(x - x)I 0 (x)]dx (3.39) 


and 


qX = -tan e sec 0e" 2lnX 


r 


A(x)e 1nX I 2 (x)dx + 2 sin 2 0I 2 (x) 


(3.40) 


where we have put 


pX 


I n (x) = 


e inX (x - x) n A(x)dx for n = 0,1 ,2, 


(3.41) 


Inserting (3.39) along with (3.29) to (3.33) into (3.21) and integrating we 
obtain 

,x 

s(l) /)- I A*(x)e~ 1nX [l 2 (x) - 2 ( x - x>I 1 <x> 


Oq 2 = tan 0 SGC 0 


- 4( x - x) 2 sin 2 0I (x) Idx . (3.42) 




Since u i y anc * W 1 are rea ^ we must set 


W <1> _ H <n * 

W 0,-2 W 0,2 


(3.43) 
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(3.44) 


and 


It Is easy to show that 




H 


(1) 

- 2,0 


<n 

2 , 


0 





cos 9 + W 


(1) 

0,0n 


sin 9 = tan 9 


Re 


pX 

A*(x)e' 1nX I 2 <x)dx 

V 

—00 


(3.45) 


( 3 . 46 ) 


(3.47) 


0 


( 1 ) 

1,1 


" 2a S l U c 


sin 9W 


( 0 ) 

lnx 


(3.48) 


and 


S,n9 -s 


2«Y 

l u c 




(1) 
1 


cos 9 + W 


(1) 

Un 


sin 9 , 


= CL 



(1) 

U 


- si wn 


( 3 . 49 ) 


where C is a constant whose specific value will not effect our analysis. 

The remaining Q's and W's can be similarly determined, but are not used in 
the following analysis. 

It now follows from (3.32) and (3.33) that (3.22) possesses a solution of 
the form 


d 2 - 


1 sin 9U' 

2 C 


n,m=-3 




<3. 50) 


where 


m * ® as n * ±0 ° (3.51) 

Substituting this along with (3.32), (3.33), (3.48), (3.38), and (3.43) to 
(3.46) into (3.22) and (3.23) yields 
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0 <2) C / H <1) _J_ s n.i/OJN 

9 1J " C V W 1J 2a S l U C W lx ; 

nJ 




0<;> * Oq'o <• tan e(w<; 


(1) u ( 1 )'\ 

♦ Vo) 


nJn 


- ^ * 


, . 2. 3 f..<0)„(l) u (0),.(l) „(0)*,,(1)1 

2 sin 9 U 0>2 - M, * K |t) V 2 _ 0 j 


2 

U S 
u cd 


OJ“ 


A U U'A- + 1A 


2a 2 sin 9 U )^ 2 A 2 


<3.52) 


where 


(1) _ (1) 
u O,2n ' y O,2 ’ 

V 0,2 n " t 0W O,2 ’ 


(3.53) 

(3.54) 


and 


(1) _ 0 (1) 
v 2.0nn ' g 2,0 • 


(3.55) 


Inserting (3.39) and (3.42), we find that 


R 0 2 E 2 sin 2 e ) = sin 6 tan fc 0( w; v/ J,- ♦ 


^ = sin 6 tan 2 e(l 


2r,f ..(0) j 2W^A __ 

1 J 3x + ^ W l n 2xx 


- W (0) J -A 
w lnn 1 xxx/ 


(3.56) 


where the J n 's depends on both x and n and are defined by 

rX 


j, n . i fit 


A*(x)(x - ?) 3 [(x - X)sin 2 e + i ^-]e" ir)X I 0 U)dx (3.57) 


J 2 = id, - J 4 ) 


(3.58) 


where 


P X 


.) 4 = 3 sin 


2 - 5 - 2 © Re 


„x 


A*(x)(x - x) 3 e" lnX I 0 (x)dx + j 


A(x)(x - x) 2 e' nX I*(x)dx 


(3.59) 
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and 


J 3 » -2 


A *<x)(x - x)|J sin 2 0(x - x)*I 0 (x) + (x - x)I 1 <x) - I 2 (x)je" 1r|X dx. 

o 

(3.60) 


Then, upon Integrating by parts and inserting (3.31), we obtain 

pX 


e“ ,n<x " x) R 0 2 d? = sin 9 tan 2 0e" iriX 


J lxx I 2 <x) " 2J 4 x I 1 (x) 


pX 


e lpX A(x)[2J 4 (x) - 2J 4 (x) + J 3 <x) - J 3 <5ojdx 


»]■ 


where 


U (1) -Q (1> 

u 2,0n g 2,0 


Inserting this along with (3.31) into 

rX 


where 


dl . 


d - i a „<o)*„(i) 

r 2,o = 2 s,n 9 aT w !n v 2,o 


and integrating by parts now shows that 

px _ ~ 

e -in(x-x) R ^ ^ 6 ~ a i s in 2 ee 1nX l*(x)U^Q(x) 


(3.61) 


Substituting (3.40) into (3.55) we find after some manipulation that 

•/(I) i _-21nx/ 9 „2inXn( 1 A 1 Q rt r A/ -,, -inx T /-. 
v 2 o “ " 2 e Va^ e U 2 oy + 2 tan 9 sec 0 [ A(x)e 1 1 ( X ) 

+ 2 Sin 2 ee" 2 lnx l 2 (x>] , (3.62) 


(3.63) 


(3.64) 


- sin 0 tan 2 0e~’ nX 


I*(x)[A(x)e lnX I 1 (x) + 2 sin 2 0I 2 (x)]dx . (3.65) 
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It is now easy to write down the solution to (3.52) by using (3.61), (3.65), 
(3.39), (3.42), and (3.47). Using the resulting formula for and 

integrating by parts shows that 


J_ 

2ir 


n ( 2 ) . 

Qi i dn 


4a 2 sin 9U 


■J G Wi * ,A ) 


X X 

+ sec 0 tan 0 J J K<x|x,x 1 >A(x)A<x 1 >A(x + ? 1 - x)dx 1 dx , 

J — 00 _oo 


where we have put 

K<x | x , x 1 ) = (x - x){[j <x - x> + si n 2 0(x - x ] )] [2< x - - (x - x) 

1 A — A -| 

It now follows from (3.8) and (3.24) to (3.28) that 


] 


- 4 sin 4 0 T ( x - x n ) 2 + (x - x)(x - x^)j|. 


(3.66) 


(3.67) 


sin 0 


n ( 2 ) , 1 , <-2.. 

°1,1 dn = 2 lS l U c 


U C A x( b 2,2 ~ b 2,2^) ^2 ( °2 A x + lA )( b 2,l 


J 2 , 1 


(3.68) 


Inserting (3.66) into this result yields 


1 A- = A + Y tan 2 © J J K(x|x,x 1 )A<x)A(x 1 )A*(x + x ] - x)dx ] dx, 

K ~oo — co 


(3.69) 


where 


1 

K 


^ u c S u c( b 2,; 



,U c U c 


(3.70) 


is the reciprocal of the scaled complex wave number in the linear amplitude 
(2.14) and we have put 


1 

Y 




i-ntr ' ' \ 
_ 

u c ' 


stu 


rc 


4a 


2 • 


(3.71) 


Equation (3.69) is the final result. It determines the amplitude of the 
instability wave. It must be solved subject to the upstream boundary condition 
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as 


(3.72) 


A -> e KX+1c P 0 


X -» -«» , 


where the real constant <p 0 is an as yet unspecified initial phase factor 
which was introduced via the arbitrary origin shifts x 0 and X 0 in (3.25) 
and (3.27). 


4. ASYMPTOTIC SOLUTION OF AMPLITUDE EQUATION 


The principal result of this paper is given by (3.69) together with the 
upstream boundary condition (3.72). The numerical solution to this problem, 
which is discussed in thenext section, appears to develop a singularity at a 
finite value of x, say _x s . In this section we determine the asymptotic form 
of the solution as x -» x s . To this end we substitute 


A = 



a 




(4.1) 


where x s and a are real constants and a is a complex constant, into the 
integral of_(3.69) and change thei ntegration variables from xi and x to 
(x s - X])/(x s - x) and (x s - x)/(x s - x) to show that 


pX pX 

K(x | x ,Xj )A(x)A(x 1 )A*(x + 5^ - iOd)^ 

J_oo '’-00 


dx = 


(X - X) 


4+1 a 


D(o) 


(4.2) 


where 


D(a) = 


3+ia 


du dv 

u 3+1a (u + v - l) 3 “ ia 


- ^(v - l) 3 + (v - l) 2 (u - v - sin 2 0[u - v + 4(u - l)sin 2 0]} 


2 _ . 2 , 




+ 2(v - l)(u - v) sin 9(1 - 2 sin"0) 


>] 


- I) 2 


E C m (v l 9) ( U - I V 

m + la V V ) 


m+i a 


m=-2 


dv 


(1 - x) 


3+ia 


^ £ m (l| l 9) \ n V - I V 

/ j m + \o v v / 
m =-2 


m+i a 


- 1 


dx , 


(4.3) 


and the C m and C m are complicated functions of the indicated arguments 
which are listed in the appendix. It is worth noting that C m is a polynomial 
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in sin 6 whereas C m is a Fourier series in 0. Carrying out the last 
integration gives 


D(o) = 


00 

-e 

n-1 


(- 1 )' 


C 1 0 (n|a) + C 1 2 (n|o)cos 20 + C ] 4 (n|a)cos 


401 


n(n + 1 >(n + 2 ) (n + 3X4 + 1 o) 


n 


C 2 q( n | a) + C 2 2 <n|o)cos 20 + C 2 4 (n|a)cos 40 
~ 1 (n + TXn + 2)(n + 3X4 + 1o> 


C 3 g(n|a) + C 3 2 <n|a)cos 20 + C 3 4 (n|a)cos 40' 
(n + 2)(n + 3)(4 + i cr) n 


(4.4) 


where (a) n denotes the generalized factorial function r(a + n)/r(a) and the 
coefficients C-| > q , . . , ^3,4 are sn ven ' n the appendix. 


Ay becomes large compared to A as x -*■ x s , and the left hand side of 
(3.69) is balanced by the integral term on the right side. Substituting (4.1) 
into the left hand side shows that the two terms will balance when a 
satisfies 


D(q) 

(3 + i a) 


tan 2 © 


1 

yic|a | 2 


(4.5) 


_Figures 1 and 2 show a and | a | respectively as functions of arg 
(l/yic), which are used to evaluate asymptotic curves In figures 3 to 6 (shown 
as dotted lines). 


Parts of the derivations in this section were done with MACSYMA (412.61 
for DEC VAX 8600 Series) batch programs. Interested readers can contact the 
second author for the programs. 


5. NUMERICAL RESULTS AND DISCUSSION 

The relevant solutions to (3.69) and (3.72) involve the two complex param- 
eters k and y, the propagation angle 0, and, the still unspecified initial 
phase factor <p 0 . But introducing the rescaled variables A/|y| ^2| K |2 anc j 
| k | x - x 0 , where x 0 and cp 0 are chosen so that 



shows that these_so1utions can be completely characterized by the two imaginary 
parameters In (k/|k|) and _ln (y/|y|), or equivalently by the arguments of k 
and y. The real part of k is the scaled growth rate of the upstream linear 
instability wave in the vicinity of its neutral stability point, and the imagi- 
nary part is the scaled deviation of the wave number from its neutral value 
corresponding to the prescribed Strouhal number deviation Sj . 

Equation (3.69) was solved numerically by using a fourth order predictor 
corrector scheme to advance the solution downstream from the prescribed 
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upstream linear state (eq. (3.72)). The double integrals were computed by 
using the Trapezoidal rule with the upstream "tails" evaluated analytically 
from the upstream linear solutions, « and y must, in general, be found 
numerically by solving the homogeneous and inhomogeneous Rayleigh's equations 
(2.19) and (2.27). Analytic solutions can only be obtained for the "tanh" 
shear layer (Huerre, 1980; 1987) where 

U c = (U (1) + U (2) )/2A, U ' - i - 1 , U£" - -2, b+ 1 r z i = 0, 

and b+ >2 " b 2 ,2 = 1 * 

Figures 3 and 5 are plots_of the instability wave growth rate | A | x / | A | 
for various_val ues arg y, arg k, and cos 9. Weon^y show results for 
-ir/2 < arg tc < 0 because (3.69) implies that A(x, <*, y*) = A*(x, k , y). 

Notice that the upstream linear growth rate is initially reduced when arg y 
is in the range -ir/2 < arg icy < ir/2. This is because the nonlinear term 
behaves like -yC 0 |A| 2 A for small |A|, where C 0 is a positive real constant. 
The effective growth rate is therefore reduced by the factor 1 - C 0 |A| 2 y, but 
(except in the special case arg y = arg « = 0) this trend is eventually 
reversed, and the growth rate rapidly increases until the amplitude becomes 
singular at some finite downstream distance for all values of y - suggesting 
an explosive growth of the instability wave there. 

This is shown somewhat better in figures 6 to 8, which are plots of the 
real part of_the scaled instability wave amplitude versus the scaled streamwise 
coordinate x. Also shown in the figures are the results computed from the 
asymptotic solution (4.1), with the singularity location x s determined from 
the numerical solution. The latter solutions clearly approach the asymptotic 
result as |xc - x becomes small. Since (4.1) implies that the asymptotic 
growth rate | A | ^ / A| behaves like (|A|/|a|)^3 i n the vicinity of the singu- 
larity, the initial scaling, i.e., growth rate = 0(e)^3, unchanged by the 
singularity. This suggests that the basic asymptotic structure of the critical 
layer will remain intact, and the present solution will not break down until 
the amplitude |A| of the external instability wave becomes order one. The flow 
will then be fully nonlinear and unsteady in the ma ; n part of the shear layer, 
i.e., it will be governed by the full Euler's equations there. 

This is quite different from the nonlinear critical layer behavior found 
by Goldstein and Leib (1988) for the two-dimensional case where the critical 
layer nonlinearity produced a redistribution of vorticity that rapidly reduced 
the instability wave growth rate to zero. The three-dimensional effects in 
the present analysis allow for vortex stretching, which completely counteracts 
this effect and produces the explosive instability wave growth. This phenomena 
is probably masked at subsonic speeds by the fact that the linear growth rate 
of the three-dimensional wave is much smaller than that of the two-dimensional 
wave. The latter wave can alter the critical-layer structure of the oblique 
waves, and the present analysis assumes that the two-dimensional wave is com- 
pletely absent in the nonlinear region. It should be possible to eliminate 
the two-dimensional wave in carefully controlled subsonic experiments, but it 
would probably be much easier to observe the phenomena at supersonic speeds 
where the most rapidly growing linear mode is oblique. While the basic ampli- 
tude equation (3.69) was derived only for incompressible flow, it applies to 
the compressible case as well, and there is no a-priori restriction on the 
Mach number. It is worth noting that the nonlinear critical layer effects 
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should be much more important at higher Mach numbers, because the linear growth 
rate rapidly decreases with increasing Mach number (Jackson and Grosch, 1988). 

Figures 6 to 8 show that the instability wave amplitude undergoes succes- 
sive oscillation upstream of the singularity. Similar behavior was observed 
in the two-dimensional analysis of Goldstein and Leib (1988) and in the calcu- 
lations of Benney and Maslowe (1975), Huerre (1977), and Miura and Sato (1978). 
The amplitude oscillations imply periodic reversal of energy transfer between 
the fluctuations and the mean flow, and possibly between the fluctuations them- 
selves. By considering the Reynolds-stress changes that occur with nutating 
elliptic vortices, Browand and Ho (1983) came up with a simple kinematic expla- 
nation for this phenomenon. The reader is referred to Ho and Huerre (1984, 
p. 410) for details. 

Viscous effects will undoubtedly alter the critical layer behavior 
described above. We decided not to include them because they would greatly 
complicate an already complex analysis. They were, however, included in the 
simpler analysis of Goldstein and Leib (1989), and we expect them to act simi- 
larly in the present situation. The former analysis suggests that viscous 
effects always delay the explosive growth and that they will eliminate it 
entirely for a certain range of values of y, k , and 9, provided an appropri- 
ate scaled viscous parameter exceeds a certain finite value. In that case, 
the solution will just go to a finite amplitude equilibrium state further 
downstream. 

While the present scaling may seem to be rather special, the composite ex- 
pansion technique of Goldstein and Leib (1988) shows that the instability wave 
adjusts to the proper scaling automatically as it propagates downstream toward 
the neutral stability point.' The nonlinear region actually sets its own loca- 
tion in the final composite expansion. The only requirement is that the insta- 
bility wave amplitude remain fairly small in the region where the linear 
growth rate is fairly small. This requirement should be very nonrestrictive 
at supersonic speeds where the linear growth rate is always small. However, 
the experiments show that nonlinearity sets in at very small amplitudes, even 
at subsonic speeds, which suggests that the nonlinear effects are localized 
and therefore confined to the critical layer (since this is the region where 
nonlinearity would first come into play; Goldstein and Leib, 1988; Goldstein 
and Hul tgren , 1 988) . 

Equations (2.7), (3.25), and (4.1) show that the explosive growth occurs 

when 


x - x s = 0(1) 


where 


x s * 


x s * x o 

1 1 /3 C 

~ e S,U 


r c 


^s in Goldstein and Leib (1988), the nonlinear solution can be consid- 
ered to be an inner solution in a composite expansion whose outer solution is 
the 1 inear wave. 
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is the singularity location in the unsealed streamwise coordinates. The 
streamwise extent of the fully nonlinear region is therefore of the same order 
as the shear layer width. 

It is also worth noting that the asymptotic instability wave amplitude is 
uniquely determined by the asymptotic solution and is therefore independent of 
the upstream conditions. Figures 9 to 1 1 show the wavelengthreduction 
y#»(Ax/A) as a function of the scaled streamwise coordinate x. The asymptotic 
results computed from (4.1) are indicated by the dashed lines. 

The authors would like to thank Dr. Lennart Hultgren for his helpful 
comments on the manuscript and Dr. S.J. Leib for helping with the computer 
program. 
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APPENDIX 


The detailed expressions for coefficients used in (4.3) and (4.4) are 


C 1 =- 


C 0 * 6 


'-1 


2 3v - 1 

v(3v - 

• 1) 

u 2 * 2(v - 1) " 

(v - 

l) 2 

7v - 

3 17v 2 

- lOv 

+ 1 

V - 

1 + (V 

- I) 2 

l 

2v - 1 

35v 2 - 

28v •» 

_5 

1 v - 1 

(v 

- I) 2 


9v - 5 

31 v 2 - 

30v ■» 

_7 

v - 1 

+ 

(v 

- I) 2 


5v - 3 

10v 2 

- 1 1 V 

+ 3 


'-2 " 2(v - 1) 


sin 0+4 


sin 0 - 4 


sin 0 + 8 


sin 0-4 


sin 0 + 4 


v - v + 


(V 

2 * 
- i r 

5v 2 

- 4v + 3 

(V 

- I) 2 

5v 2 - 

4v + 2 

(V 

- I) 2 

1 

CM 

> 

cr* 

8v + 3 

(V 

- I) 2 

3v 2 

- 3v + 1 


( v - 1)‘ 


(V - 1>‘ 


sin 0, 


4 

sin 0, 


sin 0, 


sin 0, 


(A . 1 ) 


(A. 2) 


(A. 3) 


(A. 4) 


(A. 5) 


C 2 

-10 

- 2x 

+ 4x 2 ) 

1 

2 

( 2 

- 5x 

♦ 4x 2 ) 

cos 20+| 

0- 

x + x 2 ^cos 40, 

(A. 6) 

Cl " 

-( 2 - 

3x + 

7x 2 ) + 

K 8 

- 

1 6x + 

IK 2 ) 

cos 20 - ^ 

( 4 - 

2x + 3x 2 ^cos 40, 













(A. 7) 

II 

o 

* o 

K‘- 

6x + 

19x 2 ) - 

■K 

12 

- 18x 

, llx 

2 )cos 20 + 

( 3 * 

2x 2 )cos 40, 

(A. 8) 

C -1 ■ 

. -(2- 

x + 

6x 2 ) + 

K 8 

- 

8x + 

5x 2 )cos 20 - 1(4 

+ 2x 

+ 3x 2 )cos 40, 

(A. 9) 


C -2 

■K 1 

+ 3x 2 ) 

1 - 

K 2 

- X 

+ x 2 )cos 20 + 

1 + X 

+ x 2 )cos 40, 

(A. 10) 

C l,0 


- 1o) n+2 “ 

2(- 

i a) 

n+2 + 

3(1 - 

1< ’ > n*2 ‘ 

2(2 - 

* o) n+2 * ? (3 - 

ia) n+2. 












(A. 11) 


C 2 ,0 


1 - la) 

n+2 

+ 

3(-ia) 0 - 

n+2 

3(1 - Icr) 

n+2 + 

<2 - ’ o) n+2 ’ 

(A. 12) 

C 3 ,0 

= 2(-l 

- '°>,»2 - 

7(- 

la) 

n+2 + 

13(1 
2 u 

- io) ru2 * 

6(2 

- lo) n t 2 * f< 3 - 

ia) n+2, 


(A. 13) 
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LO | OJ 


c l,2 ■ -‘- 1 - 1a> n+2 * 4< - ic,, n + 2 - 6(1 - * o) n+2 + 4(2 “ io) n + 2 ‘ <3 ' ,o) n.2- 

(A. 14) 

(-1 - i ct)„ . o - 8(-ia) n+2 + 9(1 - ia)_, - 4(2 - ia)^ ? + *<3 - ia> n+2> 


n+2 


n+2 


(A. 15) 


C, , - -2<-l - * ^<-1<7) ni , - r-d - - f<2 - '«> 


11 


'3,2 


n+2 T 2 n+2 2 


n+2 T 2 


n+2 


- 2^3 - i°^ n+ 2’ (A. 16) 

C l ,4 * I '' 1 - ' CT> n+2 - 2 <- * o ) n +2 * 3< ’ ' '°>n,2 ’ 2(2 ‘ io> n,2 * ! <3 ‘ , o ) n , 2 - 

(A. 17) 


C 2,4 = " ~ 1a) n+2 + (_1a) n+2 “ (2 " ia) n+2 + 2 (3 " 1a) n+2’ 


(A. 18) 


C 3,4 = i ( - ] * 1<7) n+2 - l ( " ia) n+2 + 2(1 " ia) n+2 “ l (2 ' ia) n+2 + I (3 ' ia) n+2> 


1 


(A. 19) 


! 
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FIGURE 1. - ASYMPTOTIC EXPONENT 0 VERSUS ARG. (1/yK) 
IN 7T RADIANS. 


27 
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